Search for: All records

Abstract Stellar streams from disrupted globular clusters are dynamically cold structures that are sensitive to perturbations from dark matter subhalos, allowing them in principle to trace the dark matter substructure in the Milky Way. We model, within the context of Λ cold dark matter, the likelihood of dark matter subhalos to produce a significant feature in a GD-1-like stream and analyze the properties of such subhalos. We generate many realizations of the subhalo population within a Milky Way mass host halo using the semianalytic codeSatGen, accounting for effects such as tidal stripping and dynamical friction. The subhalo distributions are combined with a GD-1-like stream model, and the impact of subhalos that pass close to the stream are modeled withGala. We find that subhalos with masses in the range 2 × 10⁶M_⊙–10⁸M_⊙at the time of the stream–subhalo encounter, corresponding to masses of about 2 × 10⁷M_⊙–10⁹M_⊙at the time of infall, are the likeliest to produce gaps in a GD-1-like stream. We find that gaps occur on average ∼3 times per realization of the host system. These gaps have typical widths of ∼(5–27)° and fractional underdensities of ∼(10–30)%, with larger gaps being caused by heavier subhalos. The stream–subhalo encounters responsible for these have impact parameters (0.1–1.5) kpc and relative velocities ∼(200–410) km s⁻¹. We also investigate the effects of increasing the host-halo mass on the gap properties and formation rate. 

Abstract Stellar streams in the Milky Way are promising detectors of low-mass dark matter (DM) subhalos predicted by ΛCDM. Passing subhalos induce perturbations in streams that indicate the presence of the subhalos. Understanding how known DM-dominated satellites impact streams is a crucial step toward using stream perturbations to constrain the properties of dark perturbers. Here, we cross-match a Gaia Early Data Release 3 and SEGUE member catalog of the Cetus-Palca stream (CPS) with H3 for additional radial velocity measurements and fit the orbit of the CPS using this six-dimensional (6D) data. We demonstrate for the first time that the ultra-faint dwarf Segue 2 had a recent (77 ± 5 Myr ago) close flyby (within the stream's 2σwidth) with the CPS. This interaction enables constraints on Segue 2’s mass and density profile at larger radii ( $\mathcal{O}\left(1\right)$kpc) than are probed by its stars ( $\mathcal{O}\left(10\right)$pc). While Segue 2 is not expected to strongly affect the portion of the stream covered by our 6D data, we predict that if Segue 2’s mass within ∼ 6 kpc is 5 × 10⁹M_⊙, the CPS's velocity dispersion will be ∼ 40 km s⁻¹larger atϕ₁ > 20° than atϕ₁ < 0°. If no such heating is detected, Segue 2’s mass cannot exceed 10⁹M_⊙within ∼ 6 kpc. The proper motion distribution of the CPS near the impact site is mildly sensitive to the shape of Segue 2’s density profile. This study presents a critical test for frameworks designed to constrain properties of dark subhalos from stream perturbations. 

Abstract Chemical cartography of the Galactic disk provides insights into its structure and assembly history over cosmic time. In this work, we use chemical cartography to explore chemical gradients and azimuthal substructure in the Milky Way disk with giant stars from Apache Point Observatory Galactic Evolution Experiment (APOGEE) DR17. We confirm the existence of a radial metallicity gradient in the disk of Δ[Fe/H]/ΔR∼ –0.0678 ± 0.0004 dex kpc⁻¹and a vertical metallicity gradient of Δ[Fe/H]/ΔZ∼ −0.164 ± 0.001. We find azimuthal variations (±0.1 dex) on top of the radial metallicity gradient that have been previously established with other surveys. The APOGEE giants show strong correlations with stellar age and the intensity of azimuthal variations in [Fe/H]; young populations and intermediate-aged populations both show significant deviations from the radial metallicity gradient, while older stellar populations show the largest deviations from the radial metallicity gradient. Beyond iron, we show that other elements (e.g., Mg, O) display azimuthal variations at the ±0.05 dex level across the Galactic disk. We illustrate that moving into the orbit-space could help constrain the mechanisms producing these azimuthal metallicity variations in the future. These results suggest that dynamical processes play an important role in the formation of azimuthal metallicity variations. 

ABSTRACT Carbon-enhanced metal-poor (CEMP) stars comprise almost a third of stars with [Fe/H] < −2, although their origins are still poorly understood. It is highly likely that one sub-class (CEMP-s stars) is tied to mass-transfer events in binary stars, while another sub-class (CEMP-no stars) are enriched by the nucleosynthetic yields of the first generations of stars. Previous studies of CEMP stars have primarily concentrated on the Galactic halo, but more recently they have also been detected in the thick disc and bulge components of the Milky Way. Gaia DR3 has provided an unprecedented sample of over 200 million low-resolution (R ≈ 50) spectra from the BP and RP photometers. Training on the CEMP catalogue from the SDSS/SEGUE database, we use XGBoost to identify the largest all-sky sample of CEMP candidate stars to date. In total, we find 58 872 CEMP star candidates, with an estimated contamination rate of 12 per cent. When comparing to literature high-resolution catalogues, we positively identify 60–68 per cent of the CEMP stars in the data, validating our results and indicating a high completeness rate. Our final catalogue of CEMP candidates spans from the inner to outer Milky Way, with distances as close as r ∼ 0.8 kpc from the Galactic centre, and as far as r > 30 kpc. Future higher resolution spectroscopic follow-up of these candidates will provide validations of their classification and enable investigations of the frequency of CEMP-s and CEMP-no stars throughout the Galaxy, to further constrain the nature of their progenitors. 

ABSTRACT We present a novel method for constraining the length of the Galactic bar using 6D phase-space information to directly integrate orbits. We define a pseudo-length for the Galactic bar, named RFreq, based on the maximal extent of trapped bar orbits. We find the RFreq measured from orbits is consistent with the RFreq of the assumed potential only when the length of the bar and pattern speed of said potential is similar to the model from which the initial phase-space coordinates of the orbits are derived. Therefore, one can measure the model’s or the Milky Way’s bar length from 6D phase-space coordinates by determining which assumed potential leads to a self-consistent measured RFreq. When we apply this method to ≈210 000 stars in APOGEE DR17 and Gaia eDR3 data, we find a consistent result only for potential models with a dynamical bar length of ≈3.5 kpc. We find the Milky Way’s trapped bar orbits extend out to only ≈3.5 kpc, but there is also an overdensity of stars at the end of the bar out to 4.8 kpc which could be related to an attached spiral arm. We also find that the measured orbital structure of the bar is strongly dependent on the properties of the assumed potential. 

ABSTRACT Gaia Data Release 2 revealed that the Milky Way contains significant indications of departures from equilibrium in the form of asymmetric features in the phase space density of stars in the Solar neighbourhood. One such feature is the z–vz phase spiral, interpreted as the response of the disc to the influence of a perturbation perpendicular to the disc plane, which could be external (e.g. a satellite) or internal (e.g. the bar or spiral arms). In this work, we use Gaia Data Release 3 to dissect the phase spiral by dividing the local data set into groups with similar azimuthal actions, Jϕ, and conjugate angles, θϕ, which selects stars on similar orbits and at similar orbital phases, thus having experienced similar perturbations in the past. These divisions allow us to explore areas of the Galactic disc larger than the surveyed region. The separation improves the clarity of the z–vz phase spiral and exposes changes to its morphology across the different action-angle groups. In particular, we discover a transition to two armed ‘breathing spirals’ in the inner Milky Way. We conclude that the local data contain signatures of not one, but multiple perturbations with the prospect to use their distinct properties to infer the properties of the interactions that caused them. 

Abstract The total mass of the Local Group (LG) is a fundamental quantity that enables interpreting the orbits of its constituent galaxies and placing the LG in a cosmological context. One of the few methods that allows inferring the total mass directly is the “Timing Argument,” which models the relative orbit of the Milky Way (MW) and M31 in equilibrium. The MW itself is not in equilibrium, a byproduct of its merger history and including the recent pericentric passage of the Large Magellanic Cloud (LMC), and recent work has found that the MW disk is moving with a lower bound “travel velocity” of ∼32 km s −1 with respect to the outer stellar halo. Previous Timing Argument measurements have attempted to account for this nonequilibrium state, but have been restricted to theoretical predictions for the impact of the LMC specifically. In this paper, we quantify the impact of a travel velocity on recovered LG mass estimates using several different compilations of recent kinematic measurements of M31. We find that incorporating the measured value of the travel velocity lowers the inferred LG mass by 10%–12% compared to a static MW halo. Measurements of the travel velocity with more distant tracers could yield even larger values, which would further decrease the inferred LG mass. Therefore, the newly measured travel velocity directly implies a lower LG mass than from a model with a static MW halo and must be considered in future dynamical studies of the Local Volume. 

Abstract In the coming decade, thousands of stellar streams will be observed in the halos of external galaxies. What fundamental discoveries will we make about dark matter from these streams? As a first attempt to look at these questions, we model Magellan/Megacam imaging of the Centaurus A (Cen A) disrupting dwarf companion Dwarf 3 (Dw3) and its associated stellar stream, to find out what can be learned about the Cen A dark matter halo. We develop a novel external galaxy stream-fitting technique and generate model stellar streams that reproduce the stream morphology visible in the imaging. We find that there are many viable stream models that fit the data well, with reasonable parameters, provided that Cen A has a halo mass larger than M 200 > 4.70 × 10 12 M ⊙ . There is a second stream in Cen A’s halo that is also reproduced within the context of this same dynamical model. However, stream morphology in the imaging alone does not uniquely determine the mass or mass distribution for the Cen A halo. In particular, the stream models with high likelihood show covariances between the inferred Cen A mass distribution, the inferred Dw3 progenitor mass, the Dw3 velocity, and the Dw3 line-of-sight position. We show that these degeneracies can be broken with radial-velocity measurements along the stream, and that a single radial velocity measurement puts a substantial lower limit on the halo mass. These results suggest that targeted radial-velocity measurements will be critical if we want to learn about dark matter from extragalactic stellar streams. 

Abstract Stellar streams in the Galactic halo are useful probes of the assembly of galaxies like the Milky Way. Many tidal stellar streams that have been found in recent years are accompanied by a known progenitor globular cluster or dwarf galaxy. However, the Orphan–Chenab (OC) stream is one case where a relatively narrow stream of stars has been found without a known progenitor. In an effort to find the parent of the OC stream, we use astrometry from the early third data release of ESA’s Gaia mission (Gaia EDR3) and radial velocity information from the Sloan Digital Sky Survey (SDSS)-IV Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey to find up to 13 stars that are likely members of the OC stream. We use the APOGEE survey to study the chemical nature (for up to 10 stars) of the OC stream in the α (O, Mg, Ca, Si, Ti, and S), odd- Z (Al, K, and V), Fe-peak (Fe, Ni, Mn, Co, and Cr), and neutron-capture (Ce) elemental groups. We find that the stars that make up the OC stream are not consistent with a monometallic population and have a median metallicity of −1.92 dex with a dispersion of 0.28 dex. Our results also indicate that the α elements are depleted compared to the known Milky Way populations and that its [Mg/Al] abundance ratio is not consistent with second-generation stars from globular clusters. The detailed chemical pattern of these stars, namely the [ α /Fe]–[Fe/H] plane and the metallicity distribution, indicates that the OC stream progenitor is very likely to be a dwarf spheroidal galaxy with a mass of ∼10 6 M ⊙ . 

Precise Gaia measurements of positions, parallaxes, and proper motions provide an opportunity to calculate 3D positions and 2D velocities (i.e., 5D phase-space) of Milky Way stars. Where available, spectroscopic radial velocity (RV) measurements provide full 6D phase-space information, however there are now and will remain many stars without RV measurements. Without an RV it is not possible to directly calculate 3D stellar velocities; however, one can infer 3D stellar velocities by marginalizing over the missing RV dimension. In this paper, we infer the 3D velocities of stars in the Kepler field in Cartesian Galactocentric coordinates (vx, vy, vz). We directly calculate velocities for around a quarter of all Kepler targets, using RV measurements available from the Gaia, LAMOST, and APOGEE spectroscopic surveys. Using the velocity distributions of these stars as our prior, we infer velocities for the remaining three quarters of the sample by marginalizing over the RV dimension. The median uncertainties on our inferred vx, vy, and vz velocities are around 4, 18, and 4 km/s, respectively. We provide 3D velocities for a total of 148,590 stars in the Kepler field. These 3D velocities could enable kinematic age-dating, Milky Way stellar population studies, and other scientific studies using the benchmark sample of well-studied Kepler stars. Although the methodology used here is broadly applicable to targets across the sky, our prior is specifically constructed from and for the Kepler field. Care should be taken to use a suitable prior when extending this method to other parts of the Galaxy.